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Ana tract 

Unsteady separated flow has been studied with 
the following two major objectives: (D to develop 
a methodology for active control, using forced 
unsteadiness provided by a spoiler flap, and (2) to 
characterize the development of tne separated flow, 
following the format ion end shedding of the 
starting vortex. An analysis has teen developed 
ror simulation of these flows involving temporally 
deforming boundaries. The Navier-Stokes equations 
are formulated in time-dependent general orthogonal 
coordinates which remain boundary aligned for 
arbitrary boundary deformations* The properties of 
the coordinates used are exploited to provide high 
efficiency to the numerical solution procedure. 
Results obtained for the becks tep~channel flow with 
an oscillating flap provide evidence of the 
efrectivenes3 of the flap as a control mechanism. 
The character i nation of separation is being carried 
out using a plane channel with a flap that rises 
into the oncoming boundary layer. 



introduction 

Unsteady flow separation has been the subject 
of research for several investigations* both 
experimental (e.g., Refs. 1-3) as veil as 
theoretical (e.g., Refs. M-6). The errorts are 
spurred by the interest in high-incidence 
aerodynamic flows intended to generate high lift 
but accompanied by repeated vortex formation and 
shedding. The latter phenomenon Itself is 
exploited via various control maneuvers to maximize 
lift in burst made, if not in sustained mode. 
Towards this goal, the present investigators have 
examined tne flow configuration of a backs tep 
channel with an oscillating flap [71 located in the 
primary separation zone of this vortex- dominated 
flow [63 . This thcoretloal/mimerical study is to 
be coordinated with the corresponding experimental 
investigations of Naglb and A chary a; the schedule 
of the experimental facility has not permitted that 
so far. Nevertheless, some results for this 
problem are included in the present paper, to 
demonstrate the effectiveness of the oscillating 
flap in altering the instantaneous separated- flow 
features of this configuration. 

A focus of the present paper is to analyze the 
process of separation itself , so that the flow- 
control meohanlen may be triggered appropriately 
and vehicle control maintained* For this purpose, 



the simpler configuration of a plane dmnnel v*l;i a 
moveable flap on its lower wall has b*ei examine! . 
For small flap-height to channel-width -stio, this 
configuration approaches that of a flap on a flaL 
plate - a configuration for which det&ued 
experimental data is expected to be available 
shortly from Acharya C93- 

The problem is formulated using the full t '.lie- 
dependent Nevier-Stoices equations in terms of tils 
vorticity « and the stream function 9. Direct 
numerical simulation of the unsteady flow is 
achieved by "direct solution of the a tream-f unct I on 
equation using block-Qaussian elimination (BGE) for 
every time-step advance of the vorticity field 
using the alternating-direction implicit (ADD 
method for the vorticity e Qua t ion HO.U The 
coordinate system employed consists or general 
orthogonal coordinates generated by 2 ggneralltid 
conformal mapping C11], with appropriate 
contraction mappings. The latter ser^e to map the 
doubly infinite channel in the physical plane t: a 
bounded region in the computational plane, as will 
as to appropriately place the uniformly spaced 
computational -grid points so as to provide 
resolution of the various length scales in the 
physical plane. The overall coordinate 
transformation is capable of generating bounder j - 
aligned coordinates for arbitrary or 1 station 01 
the flap. To accommodate a moving flap, the 
coordinate transformation necessarily varies wim 
time. Consequently* the transformatlor metrics In 
the flow equations are tlme^varying. The edver: e 
effect of this on the direct solution of the sti earn 
function equation can be significant. Rowever, as 
described in Refs- [73 and Ci2], the situation can 
&e very effeotively treated by exploiting the 
properties of conforms 1 mappings, thereby 
maintaining the efficiency of the adi-ugE nuroer cal 
technique. 

The physics Of vortex-dominated unsteady f ow 
displays widely varying length and tiin<» soaleB. 
Resolution of these very disparate scales requii-es 
an enormously increasing number of grid points .is 
the Reynolds number Re increases, with extremel;- 
lerge computational times. Therefore, direct 
simulation of such flows for higher ft» is facet 
with computer storage limitations for obtaining 
accurate resolution of the vorticity f;;eld wnld 
directly affects the pressure - a paraneter 
conveniently monitored experimentally. The us<> of 
a turbulence model may relieve some of these 
limitations. But the ambiguity associated wiu 
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turoulenoe modelling raises some concern p specially 
when characterization of the instantaneous flow 
field is of interest. It is recognized, however, 
that this congern is greatest for steady 
configurations exhibiting unsteady characteristics, 
i^. , ror flows with self-induced unsteadiness, 
for flows with forced unsteadiness, soon as the 
configurations with a moving flap, the use of a 
turbulence model can yield meaningful unsteady- flow 
results, on the tine acale of the forced 
unsteadiness. 



Mathematical Formulation 

governing Differential equations 

The time-dependent Navier-Stokes equations 
expressed in terms of the stream function ^ and 
vorticity us with variable effective viscosity v e , 

and time-varying general orthogonal coordinates, 
can be derived as follows. 

Vortioity-Transport Equation 

In general vector form and physical time t, 
tne transport equation for vorticity u is given as 

|2 - «-<W) * V'<VJ) = VxCV-<v e der v)] (i) 

where Z - vxV > (2) 



V-Cv VV) - .W - v^Vxw , 
o se 



(7) 



so that, using Eq. (3). the curl of tna second t«:-m 
on the RHS of Eq. (7) leads to 



« h jh 2 3^ n, 



For the first term on the HHS of Eq, if is 

noted that, for 2-D flow. 
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For 2-D flow in the plane of the unit vectors (e^ . 

e 2 J perpendicular to the unit vector Zy is the 

only non-zero component ot u and ts denoted as <n. 
For suah flow, y and def V can be written as 
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Consequently, 
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where 
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In view of Eq. C4). the RMS of Eq. (1) consists of 
two major terms: the second or these is considered 
first, Voting that 
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its curl becomes 
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To deal with the regaining terms contained 
within the RHS of Eq. (I), it is noted that 



Forming the curl of Eq. (11) c oapl«t*6 th«? 
derivation or the required expression rcr the nH : 
or Eq, CD, which then hae the following rinai r:rm 
in the computational coordinates (t^C^Oi 
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w!«ire A 1# $i 2 are *a given by Bqs. (13) and (13)* 
ai* g Q , t; ft represent grid-velocity components ♦ 



At stationary non-porous walla , the boundary 
values for * and u are obtainable from Eqs, (19) 
and (20)* with v being set to zero* 
2 

The Inflow and outflow boundaries are placed 
at — 1 and respectively. The esymptotia 

form of the governing equations (Hi) and (17) » as 
are solved numerically to determine the 

values « w (Cg) and *_(5 2 > at these boundaries. 

Initial Conditions 

At the boundaries* the initial conditions are 
exactly the boundary values described in the 
preceding section. In the interior, the initial 
conditions are set; according to the following 
relations i 



.Equation (iH) is the required vortieity~ 
ti insport equation for variable-viscosity flow. It 
siouid be noted that* if the fluid viscosity were 
c< istant« v.hc last term In the RHS of Eq. (1M) 
Vi.iishas, while the remaining two terns simplify to 
y Bid the Ocrrect expression far constant-viscosity 
f: >4. 

Si. «aam-Pun<stlon Equation 

This :la obtained by introducing the definition 
o:' the atroaa function, namely. 



h 1 *i 



'1 1« 2 3^ * '2 
into the definition of vortieity given as 
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<16) 



n^aae, the stream function equation is obtained as 
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and 
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Thus, the initial distribution for * consists 
of the inflow/outflow boundary values of the stream 
function, while that for u corresponds to the 
asymptotic values scaled by the local metrics. 



Outline of Solution Procedure 

Following the initial setting up of the 
mathematical problem, the general solution 
procedure consists of throe major steps; 

& Or id generation 

e Solution of stream function equation 

© Solution of ^ortioity transport equation. 

These steps are described here briefly; additional 
details are contained in Refa. [73 and Cl2]. 



Bo undary Conditions 

On tits moving flap , the boundary velocity is 
mn-zero, l,e M 



wall 



Vl * V,V 



Uie of Sqs. <15) and (16) leads to the following 
conditions for y and ui, respectively, at the moving 
flap: 



Wi 
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wall. 
(20) 



At each disorete tine instant of the unsteady- 
flow calculation, a boundary aligned coordinate 
system la generated using a generalized Schwarz- 
Christoffei conforms! mapping technique. These 
coordinates are appropriately distributed using 
suitable one-dlmenslonal clustering transformations 
so ss to map the doubly infinite channel 
configuration to a finite computational region, 
while simultaneously attempting to provide 
resolution of the various Important length scales 
of the flow. The oonformal mapping varies with 
time, so as to accommodate the moving flap; the 
clustering transformations are maintained 
independent of tiae, so as to retain the efficiency 
of the direct-solution procedure used for the 
stream function equation. 

The stream function equation (17) ie a Poisson 
equation and Is rearrange4 such that the 
differential operator which operates on * is time 
independent, and all the time variation is 
contained totally within the source term of this 
Poisson equation. The equation is disoretized 
using central differences and the resulting system 
of linear coupled algebraic equations is solved 
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using the direct blocK-Gaussian elimination (BCE) 
method developed by Osswald et al. [10], 

The temporally parabolic, spatially elliptic 
vortioity-traneport equation la solved very 
efficiently by an alternating-direction implicit 
(ADI) technique r with implicit treatment of the 
wall boundary condition for vortioity. For those 
cases planned to be tested using turbulence 
modelling, the Baldwln-Loroax [133 turbulence model 
la to be employed, it is recognized that the model 
haa deficiencies for the unsteady separated flows 
under consideration- With careful interpretation 
of the model r the present investigators have been 
able to treat small separated regions [ml. In a 
concurrent effort , various algebraic models as well 
as differential-equation models are being examined 
so as to determine if any of these other models 
perform superior to the Baldwln-Lomex model. 



Results and Discussion 

Bgckstep Channel with Control _R_ap 

This configuration (Fig. 1) enables the 
examination of flow control induced by a surface- 
mounted flap introducing spanwise vortioity into 
the separated- flow region. For sufficiently low 
Reynolds number, a steady-flow separation is 
Observed, and the reattachment length can be 
significantly reduced via the flap mechanism. As 
the Reynolds number is incr eased , spanwise 
vorticlty connected Into the separated- flow region 
13 not all connected out Of it. This leads to 
accumulation of vortioity, resulting in a large, 
spanwise vortex which is eventually shed, carrying 
the accumulated vorticlty with it downstream]. The 
process repeats, and the flow is pereistently 
unsteady, the unsteadiness being enhanced byj 
sympathetic vortex formation and shedding along the 
upper wall as well* Controlled shedding of jthe 
large vortex requires controlled removal of ; 
spanwise vorticlty from it so aa to control p,he 
vorticlty accumulation inside It. An oscillating 
flap can serve to achieve this as well. The) 
configuraiton was investigated experimentally by 
Reisenthal [153. The Reynolds number values! used 
in the experiment were rather high to enable] 
accurate spatial resolution in the direct numerical 
simulation. Values of He up to 32,000 (bases' on 
channel width) have been considered in the present 
work , 

Figure 2 shows the instantaneous stream- 
function contours for the backstop channel, of unit 
height, with a rlap of height 0.25, located at' x=»1 
and oscillating einuaoi dally "at reduced frequency 
F=0.l6 based on channel height. This frequency 
corresponds to a reduced frequency of O.Qfc based on 
flap height. The Reynolds number Re, based on 
channel height, Is 2000. Comparison of these 
results with those obtained previously Osswald et 
al. [16]) for the pacVstep channel without the flap 
shows that, for example, at Time * 60, the control 
flap leads to clear reduction in the reattachment 
length of the primary separation region. 

Associated with the control flap are various 
parameters such as Its location and geometry as 
well as the waveform and frequency or Its motion. 
Of these, the frequency is considered to be the 
more significant one* Figure 3 shows the results 
obtained for the case with F~0.25, i>e., a reduced 



frequency of 0.0625 based on flap heigh--. The 
length of the primary separation zone is |geperall| 
smaller than for the case with F-0.16 in Fig* 2. 

In order to characterize the formation and 
shedding of the starting vortex and Lhe iieparaticii 
that follows, the configuration considered is the', 
of a flap that rises up to a designated ; ingle lnt-> 
a boundary layer and then remains statiojiary unti 
the flow stabilizes. Figures *l and 5 show the 
reaults for such a flap on the lower wall of a 
plane channel for Re=»2000 and BOO0, respectively. 
These values correspond to Re f -500 and 2000. 

respectively, based on flap length. T*e case of 
Re^OOO is shown in Fig. 6. Also, for this las': 

case, the channel width is 8 flap heights, so thsit 
Re -32, 000, and is intended to better represent ths 
flat-plate configure iton used in the oor respond it: 5 
experiments of Koga [17] and Acharya [«j! 

It is necessary to mention, however, that t"s 
experiments of Koga [17] are Tor Re f *500*000. T.;e 

ll 

computational results described in the present 
paper have been obtained using direct nijtuiation of 
the unsteady Havi er-S tokes equations, l-e>, with no 
turbulence modeling. As many as (598x51 ) grid 
points have been employed in these commutations. 
Even with such refined grids, with the best 
clustering possible, the direat solutions have bi en 
limited to Ro^H.OOO- On the other hancV, the 

lowest value of Re r at which reliable experiment! 1 

measurements may be possible has been estimated is 
10,000 by Acharya [93- Hence, the numerical 
analysis and computer code are present L>| being 
generalized to accommodate a turbulence model, s<- 
that higher values of Re f may be considered in zi 0 

computations. At this time, the turbuL<^c& mode, 
of Baldwin and Lomax [133 is planned to be used. 
Reviewing Eq. <1U), it is observed that the last 
term on its RRS may not be expressed totally In 
terms of vorticlty alone. The appearai;e of 
velooity derivatives in this term furt'tor 
reinforces the 0<At) accuracy of the p^sent 
calculations. Therefore, effort needs to be 
directed to improving the temporal accireey of the 
numerical scheme as has been done in fisl\ [183- 



Concluolon 

A capability has been developed for 
numerically simulating 2~D flows in tenporally 
deforming geometries » 7he procedure depends, ir an 
essential manner, on recognizing that, with tlmr* 
independent grid clustering, the stream function 
equation can be rearranged to have a tluc- 
Independent differential operator, even in the 
presence of time-varying boundaries. The analysis 
has been demonstrated via preliminary results 
obtained for the backstop channel with an 
oscillating flap as an active flow-control 
mechanism* The frequency of the flap motion is in 
important parameter in this mechanism. 

Further analysis of the separation process 
itself is being examined via a flap or: a plane- 
channel wall, approximating flow over a riat pi.;te 
with a flap. A turbulence model is being 
incorporated into the analysis so as to obtain 
solutions for the higher Reynolds numbers u3ed in 
the experiments. 



Separated flows with forced unsteadiness are 
far *or© reproducible then those with self- induced 
unsteadiness as the latter display chaotic-flow 
characteristics. The aerodynamic rorces associated 
with unsteady separated flows have peak values 
many-fold larger, though transient, than those In 
static configurations. Understanding and control 
of such flows can play a vital role in achieving 
euen aerodynamic enhancements as needed in 
aupermeneuverabillty of aircraft. 
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Fig. 2. Unsteady Flow in Back-Step Channel with an Oscillating Flap 
Re = 2000, F * .16, 45° < 8 < 90°. (375,34) 
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Fig* 3. Unsteady Flow in Back-Step Channel with an Oscillating flip 
Re 2000, F= .25, 45° < 6 < 90*, (375.34) 
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